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I. INTRODUCTION 

O 

\ According to the Hohenberg-Kohn theorem [1[ , the ground state of an interacting electron 

gas in an external potential is a unique functional of the density. Together with a practical 
method of prescribing the density [2(], a density functional theory (DFT) has played a vital 
role for calculating the ground state energy of electron systems. Whether or not the DFT 
can be applied to a nucleus which is a self-bound system is an important question. Since 
the nuclear Hamiltonian includes no single-particle external potential, it is obvious that the 
proof of the Hohenberg-Kohn theorem does not apply for the nuclear ground state even 
though the application of the DFT is justifiable. 

There are several papers which appear to support the DFT for nuclei jl, 0, [1]. The 
^ ' arguments made in these papers all assume some sort of intrinsic density for which the DFT is 
discussed. For example, the intrinsic density is formed by putting the center of mass motion 
^ , in some potential well or by assuming a symmetry violating intrinsic state. In the former 
^ \ case one has to separate the genuine internal motion from the center of mass motion, which 
is in general not trivial as in the case of large space shell model calculation encompassing 
major shell mixing. In the latter case a physical wave function is obtained by restoring the 
symmetry by a projection procedure. This approach has however only a limited validity, 
that is, assuming an intrinsic state is already an approximation to a many-body theory. An 
intrinsic shape, if it is meaningful at all, should appear automatically from a theory which 
has no recourse to the existence of such an intrinsic shape jl, 0] ■ The two-a cluster structure 
for *Be just comes out from a calculation which involves no such assumption f\. Our recent 
four-nucleon calculation [8i] has succeeded, without assuming a cluster ansatz, to show that 
some of the excited states in ''He have 3iV + N cluster configuration. 

A generalization of the DFT is discussed by introducing a pair density [ol, [13] as a key 
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quantity to characterize the system of interacting many particles. The pair density or two- 
particle density gives a deeper insight into the internal structure of the system, especially 
into the correlated motion. The purpose of this paper is to examine internucleon correlation 
functions (iCF) since the energy of the nuclear ground state is manifestly a functional of these 
functions. Following the Hohenberg-Kohn theorem, we can unambiguously prove that the 
nuclear interaction can be uniquely determined by the iCF, that is, the nuclear Hamiltonian 
and the iCF has a one-to-one correspondence. Examples of iCF are given for s-shell nuclei. 
They are calculated using accurate wave functions obtained with realistic interactions. 

Since the nucleon-nucleon interaction depends on the spins and isospins of the nucleons, 
we have to consider the iCF in different spin-isospin channels. We discuss the relation 
between the various terms of the nucleon-nucleon interaction and the iCF as well as the 
asymptotic forms of the iCF. Information on the iCF is expected to be important for a class 
of variational calculations which use correlated trial wave functions including correlation 
factors such as variational Monte Carlo llj, coupled-cluster theory [l2|, Fermi hypernetted 



chain theory [13|, and cluster expansion method [Ij] and for a many-body theory using a 
unitary transformation of the nucleon-nucleon interaction (l5| . 



II. ENERGY AS A FUNCTIONAL OF INTERNUCLEON CORRELATION 

FUNCTIONS 

A. Definition of internucleon correlation functions 

A Hamiltonian for a nucleus with nucleons is taken as 

^ 1 1 

H = K + V = y — p2 - ^ y- /^x 

i=l i<j 

where m is the nucleon mass, P = ^^Pi is the total momentum, and the center of mass 
kinetic energy is subtracted so as to calculate the internal energy of the nucleus. 
A nucleon-nucleon interaction Vij may be expressed as follows 



(2) 



p 



where rij = \rij\ with = ri — Vj being the relative distance of nucleons i and j. Three- 
body forces are ignored for the sake of simplicity. The operators O^^^ denote various terms 
of the nucleon-nucleon potential. For the first eight terms, e.g., they are defined as 

Off = S,,, Off = S,,T,-T„ Off = iL.S),,, off = {L.S),,T,-T,, (3) 

where Sij = 3(fy ■ cri){fij ■ (jj) — CTi ■ a-j is the tensor operator, and {L ■ S)ij is the spin-orbit 
operator where L = x iVij with iZij = —i[d /drij) and S = \{(Ti + cTj). The Coulomb 

potential is included in Eq. ([2]) with /vij and off = PinFj-K where Pjjr is 1 for 

protons and for neutrons. With the use of Pi — pj = 2/j,7rjj, the internal kinetic energy is 
rewritten as 

Kj Kj 
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Let \1/ and E denote the wave ftmction and energy of the ground state of the nucleus, that 
is, H"^ = E"^. The wave function ^ satisfies all of symmetry properties such as translation- 
invariance and rotation-invariance. Assuming that ^' is normalized, the energy is given as 
the expectation value of H, E = {^\H\^). Because ^ is antisymmetric with respect to an 
exchange of nucleons, E is reduced to 

N{N - 1) 



E^ 



(5) 



The expectation values of the kinetic and potential energy terms can be expressed in terms 
of the iCF. The kinetic energy term reads 

1 77^2 1*) = j dkk^^\5{ni2 - fe)|*) = dkk^C{k), (6) 

where k = {k, k) is just an integration variable, and C{k) is the iCF in a momentum space, 
which is defined by 

c(k) = m^p^m- (7) 

The calculation of C{k) is easily performed if the wave function ^ is given in the momentum 
space because tvu = {tti2, t^u) is then just a multiplying operator. 

In exactly the same way, we can express the potential energy term as 



{^i/\vi2\^) = ^ / drv^\r){^\5{ru-r)0^^\^) 

= V / drr\^\r)C^\r), (8) 

where C'^^r) is the iCF corresponding to the operator of type p in the nucleon-nucleon 
potential 

C'W(^)^(^|^(!H_ZI)oCp)|^). (9) 



Here r = (r, f) is not a dynamical coordinate but an integration variable. The iCF for the 
Coulomb potential is defined similarly 

C^coui(r) = (^| ^(^^^-^) p,A|^). (10) 
The energy E is thus manifestly a functional of several scalar iCF, C{k) and C^'P\r): 

E = {N-l)^J^^dkk^C{k) + ^^°°drrV^)(r)C(^)(r). (11) 

As seen above, the energy can be expressed in terms of iCF. They are different from a 
two-nucleon density, which is defined as 

pir,R) = {^\S{ri-XN-R-lr)S{r2-XN-R+^r)\^) 

= (*|5(ri -r2- r)S{l{r, + rs) - xn - H)|*), (12) 
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where a^Tv is the center of mass coordinate of the nucleus. Because vu is independent of 
where the center of mass of the two nucleons relative to the total center of mass is located, 
i.e., of the coordinate |(ri + — x^, in calculating the expectation value of the potential 
energy we can integrate over R, that is, we only need 

J dRp{r,R), (13) 

which is nothing but the iCF, C^^\r). Note that C^^\r) is different from an intrinsic 
one-body density, which is defined as 

p'-^^r) = N{^\6{ri-XN -r)\^). (14) 



B. Extension of Hohenberg-Kohn theorem 

We think that no consensus has yet been reached on the existence of DFT for the nuclear 
Hamiltonian ([T]). We can prove, however, that the iCF can constitute a set of basic variables 
for the nuclear system. Obviously C^^\r) for the ground state are uniquely determined by 
Vij , and hence they are functionals of Vij . Following the proof used in Ref . [jj , we can prove 
that, conversely, Vij is a unique functional of C^^\r). For this purpose we only need to 
show that Vij is uniquely determined by C^^^r). Let us assume that the ground state of the 
Hamiltonian ([1]) is non-degenerate. Assume that, contrary to the statement to be proved, 
there is another potential v'^j, which gives rise to a ground state wave function \1'' and an 
energy E', resulting from the same iCF C^^\r). Clearly "if' cannot be equal to \E', because 
they satisfy different Schrodinger equations. Let H' denote the Hamiltonian obtained by 
replacing Vij with v'^j. Then, from the Ritz theorem, we have that 

E' = < = {^\H + V' - (15) 

Here the inequality (^'|if'|\l'') < (\l'|if'|\l') holds because is different from Using 
Eq. ([n]) leads to 

E' < ^(^- 1) r drr^[v'^P\r)-v^P\r)]C'-P\r). (16) 
2 p Jo 

Interchanging primed and unprimed quantities yields the result 

E < r drr^[v^'p\r)-v'^'P\r)]C^^\r). (17) 

2 p Jo 

Adding up Eqs. (fT6l) and (ITTl) leads to the well-known inconsistency 

E + E' <E + E'. (18) 

Thus we can conclude that Vij is a unique functional of C'^^\r). Since Vij specifies H 
unambiguously, it is concluded that the nuclear ground state is a unique functional of C'^^\r). 
The ground state energy E takes a minimum for the exact iCF. 

The wave function ^> depends on 3A^ — 3 variables as well as the spin and isospin coor- 
dinates. It is therefore hopeless to obtain an accurate wave function for > 10 using a 
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basis expansion method. Contrary to this approach, the above consideration tells us that to 
calculate the ground state energy accurately we need to know about 10-20 iCF which are all 
single- variable scalar functions. It is interesting to know the characteristic behaviors, e.g., 
the shapes and magnitudes of these CF. 
The iCF satisfy the following equations 



oo 

2/ 



dkeC{k) = 1, 

poo cy 

I drr'C^'\r) = ^^^_^^ {^\J20^^\^). (19) 



Using the identity 



N 



Y,{n-r,f = NY^{r,-XA)\ (20) 

i<j i=l 

the root mean square matter radius of the nucleus can be calculated from a moment of 

C(i)(r) as 

{^.Y^ir^-x^f) = ^ rdrr^C(^\r). (21) 
i=i 

An interesting relation is obtained by expressing the left side of Eq. ( l2Ti) with the use of 
the one-particle density 

dr r^p^^^ (r) = j dr r^C'^^^ (r ) . (22) 

For a spherical density we have the following relation 

/•oo AT _ 1 poo 

/ drr^p'^^^r) = - drr^C^^\r) (23) 

with 

p«(r) = j drp^'\r) = iV(vI/|^ltiZ_^^^d^|vI/). (24) 



C. Internucleon correlation functions in spin and isospin channels 

The characteristics of nucleon-nucleon potentials may be more transparent if we decom- 
pose them into four spin and isospin channels of two nucleons, (ST) = (10), (01), (11), (00), 
instead of using the operator representation of Eqs. ([2]) and ([3]). To do this we use the 
following identities 

1 = + + ^ + pp, ^1-^2= p^ - 3/f ^) + - 3/f °\ 

ri-r, = -3ifnifni^r^-3if^ 

ar-cT,Tr-T, = -SP^ - SP^ + P^ + gP^ , (25) 
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where Pif^^ is the projection operator which projects onto the state with {ST) of the two 
nucleons i and j. The nucleon-nucleon potential ([2]) is expressed as 



v^^^ir,,)0,,Pp, (26) 

(ST)0 



where the summation label O indicates summing over various components of the nucleon- 
nucleon interaction. For example, they stand for central (1), tensor (T), and spin-orbit 
(LS), and the corresponding operators Oij denote 1, Sij, and {L ■ S)ij, respectively. For the 
potential ([2]), the form factor v^'^^ is related to those of f as follows: 

^;(ll)=^;(l)+^;(2)+^(3)+^(4)^ ^H) = ^(5) + ^(6) ^ ^(H) ^ ^(7) ^ ^(8) ^ 

^;fo)=^;a)-3^;(^)-3^;(=^) + 9^;W, v^^'^ = v^'^ - 3v'-'\ v^^'> = v^'^ - 3v^'\ (27) 
The iCF terms corresponding to the potential form of Eq. (1261) are 

Cr\r) = m^-^^^^0,,Pir'm. (28) 
The relationship between the two iCF, Eqs. ([9]) and ( l28l) . reads 

Cfo) = l(3C(i) + - 3C(=^) - C(^)), 4'°) = ^(C(^) - C(6)), 
4f = 1(C(^)-C(«)), 

cf) = l(3C(i)-3C(^) + C(3)-C(^)), Cf)=0, C£^^ = 0, 

= l(9C(i) + 3C(2) + 3C(3) + C(^)), = ^(3C(5) + C(6)), 

16 4 

43^) = 1(3C(^)+C(«)), 



LS 
(00) 



_L(c«-C(2)-C(3) + C(^)), 4°) = 0, 4s°^ = 0. (29) 



Equation (jTTj) is rewritten using these zCF as follows 

-2 /-oo 



= (iV-1)- r dkk'C{k) + ^^^~^^ V r drr^v^^^\T)&^^\r). (30) 
Jo 2 ^^^^^ Jo 



III. SPECIFIC EXAMPLES 



Recently we have developed a method of calculating matrix elements for the interaction 
of Eq. ([3]) as well as various types of iCF using correlated Gaussian functions with the 
orbital motion being described in two global vectors 1^. The accuracy of the formulation 



has been tested by comparing to other calculations for iV = 3 — 4 nuclei. An application of 
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the method to studying excited states of ^He has met a fair success, reveahng an inversion 
doublet picture arising from ^H(t)+p and ^He(/i)+n cluster structure [st]. 

The wave functions are expressed as a combination of many basis states, each of which 
has the following LS coupling form 

^iLS)jM,TM. = A[^pt^-vr%MATr^ , (31) 

where the square bracket [. . .] stands for the angular momentum coupling, and the antisym- 
metry of nucleons is met by the antisymmetrizer A. The spin and isospin parts are expanded 
using the basis of successive coupling, e.g., 

i^SM^ = |[' ■ ■ [[[illsiailsias] ' ' 'l^Ms), (32) 

where the set of intermediate spins (5*12, >S'i23, . . .) is allowed to take all possible values for a 
given S. 

The orbital part ipi^f^^^ is given as follows 

F^LiL2)lm{ui,U2,A,x) = exp (^-^xAx^ [yLA^ix)yL2{u2x)]LM, (33) 

with yLAiir) = r^YLM{r), where x = {xi, X2, . . . , a^Ar-i) is a set of relative coordinates, say 
the Jacobi coordinate set, and Ui and U2 are (A^ — l)-dimensional column vectors which define 
the global vectors. Here xAx = ^lAijXi ■ Xj with Aij = Aji and uix = ^iJi^ uuXi. 
As we see, each basis function is characterized by a set of parameters, A = (Aij), ui, U2, 
Li, L2, L, S12, S123, ■ ■ ■ , S, T12, T123, . . .. 

The calculation of Hamiltonian matrix elements is made possible with the aid of the 
generating function g 

g{s; A, x) = exp ^ — -xAx + sx^ , (34) 

where s is an (A^ — l)-dimensional column vector whose ith element is a 3-dimensional vector 
Si. By expressing Sj with 3-dimensional unit vectors ei and 62 as Sj = XiBiUi. + A2e2U2i, 
the basis function fl33|) is generated as follows: 

F{LiL2)lm{ui,U2,A,x) = ^jj^ J J deide2 [YL,{ei)YL2{e2)]LM 

X Q^L,Q^L2 9{^ieiUi + X2e2U2;A,x) , (35) 



Ai=A2=0 



where 



Formulas for the matrix elements are given in Ref . 1^ . In Appendix [X] we giv e a formula 
to calculate the zCF for the spin-orbit force, which was not included in Ref. [161 ]. 

We study the iCF of s-shell nuclei, d, t, h, and "^He (a). We also show the iCF calculated 
for the first excited O"*" state of ^He, which is called a* in this paper. Although this state 
decays into the t + p channel with a width of 0.50 MeV, approximating it as a bound state 
is fairly good, and the iCF of a* are calculated using the wave function obtained in that 
approximation jl]. Because a* is a spatially extended state with a cluster structure of t + p 
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and h + n, comparing the iCF between a* and other cases reveals how much the shapes and 
magnitudes of the iCF are modified by the structures of the underlying nuclear states. 

We use the AV8' potential jl7|] and the GSRS potential [ll] as the two-nucleon interaction. 
Both of them contain central, tensor and spin-orbit terms. The L"^ and {L ■ Sy terms of the 
GSRS potential are ignored. The radial form factors v^^'^^ of the two potentials are displayed 
in Fig. [TJ As is well-known, the longest-range attraction of the two-nucleon interaction is 
that belonging to the term. In the intermediate region (1 < r < 1.6) fm, the singlet 
even central potential is most attractive, and then the v!^^^ and vf^^^ terms follow. 
The central potentials all have strong short-range repulsion. Generally speaking, the GSRS 
potential is softer than the AV8' potential. The tensor force of the latter potential is much 
stronger at r < 1 fm than that of the GSRS potential. The GSRS potential has no f Lg°^ 
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component. In spite of these differences, the two potentials give rather similar binding 



energies for the s-shell nuclei [16 . 

Figures [2] and [3] compare the iCF, Cq 



[r), calculated using the wave functions obtained 
with the AV8' Hamiltonian with those obtained with the GSRS Hamiltonian. The amplitude 
of c[^^^ is by far smaller than the others because of the purely repulsive nature of the 
corresponding potential vf^^\ and thus it is not shown in the figures. Common to the 
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FIG. 2: The iCF for a, a* (the first excited 0+ state of ''He), t, h, d calculated using the AV8' 
potential. Note that different scales are used for the vertical axes. 
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two cases is that the three iCF, Cr^ , C\ and C[ , among others have much larger 
amphtudes than the others, all having a peak at around 1 fm. This is understandable from 
the characteristics of the Vq'^^ curves shown in Fig. [H The more attractive the potential 
component, the larger the amplitude of the corresponding iCF. The central potential for 
even partial waves has a minimum around 1 fm. Since the central and tensor forces in the 
triplet even channels couple, the c!^^^ curve also has a peak at almost the same position 
as that of C[^^\ though the v!^^^ potential does not have a minimum around Ifm. Both of 
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FIG. 3: The same as Fig. [2] but calculated using the GSRS potential. 
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C[^^^ and d^^'^ curves have a peak around 1.4-1.5 fm because the central and tensor forces 
in the (ST) = (11) channel couple and the central force has a shallow attraction beyond 
Ifm. The peak position of the iCF for the Coulomb potential is the same as that of the 
main terms of the triplet even channel. 

(ST) (ST) 

Because of the strong short-range repulsion of the central potentials , all the Cq 
vanish near the origin. The spin-orbit force (t'Ls^'*) has a very strong attraction around Ifm 
but it is confined to the short distance region. Comparing Figs. |2] and [3l the difference of 
the iCF produced by the potential models is mild in spite of apparently different behaviors 
of the potential form factors, for example of the tensor terms. 

The s-shell nuclei all produce similar shapes for each of the iCF. It is interesting to 
note the iCF of a* exhibit the patterns similar to the other cases, despite the fact that its 
structure is quite different from that of a. Of course the amplitudes at larger values of r are 
much larger for a* than for the other cases. 

A remarkable characteristics of the iCF is that their asymptotics for a given nucleus are 
the same for all the Cq'^\ For example, the eight iCF calculated for a using the AV8' 
potential follow ~ exp(— 2Kr)/r^ °^ with k = 0.75 fm~^ for r larger than 6fm. 

The case of the deuteron is easily understood. The deuteron wave function consists of 
the S- and D-wave components 

^ = (t)Q{r)YQQ{r)xiMrim + </)2(r)[y2(T^)xi]iA/?7oo, (37) 

where xi ^i^nd r/o are the spin and isospin functions of the deuteron. For large values of 
r for which the nuclear potential between the two nucleons in the deuteron are negligible, 
the Hamiltonian for the deuteron reduces to the kinetic energy alone, and thus the radial 
function (/)^(r) should be given by a solution of the free-particle Schrodinger equation with 
the negative energy {—h'^K^/m) of the deuteron, that is the spherical Hankel function of the 
first kind hf\iKr). The asymptotic form of the deuteron wave function is therefore given 
by 

^=J2 Kihf\iKr)[Ye{r)xi]iMmo (38) 
e=o,2 

with suitable coefficients Kg. The iCF of the deuteron for large r reduces to 

Co^'V)- Yl K,Kp{[Y,ai]iMVoo\OP^'''^mxi]iMVoo)h'i\tf^r^^^^^ (39) 

£,^'=0,2 

All of the zCf for the deuteron displayed in Figs. [2] and [3] satisfy the above behavior for 
r > 5fm. 

The asymptotic behavior for other cases is discussed in Appendix[B]by taking into account 
the Coulomb force. The case of a is understood by taking the nucleus R (in the notation 
of Appendix [B]) as the pn system, which gives Z = 2. See Eq. (lB5p . As discussed in 
Appendix [B| the asymptotic behavior is given by r^^^^^ exp(— 2Kr), where r] defined in 
Eq. (lB7p becomes 0.046 for k, = 0.75 fm~^. Thus we can understand the asymptotic behavior 
of the iCF noted above. 

The iCF C{k) for the kinetic energy operator is shown in Fig. |H The AV8' and GSRS 
potentials give qualitatively very similar results. The behavior of C{k) for small values of 
k is given analytically as explained in Appendix [Bl The numerical results confirm that 
the asymptotic form, Eq. flB12l) . agrees with C{k) for small values of k. The behavior 
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FIG. 4: The iCF for the kinetic energy operator. 



of C{k) for large values of k primarily reflects the short-range central correlation and the 
tensor correlation involved in the wave function. The enhancement of the curve around 



k ~ 1.3 fm~^ is due to the tensor force, as shown in Ref. 191 . 

Figure displays the kinetic energy density of the two-nucleon relative motion, 
{h'^/m)k^C{k), as a function of k. We clearly see that the tensor correlation increases the 
kinetic energy density beyond k = lfm~^. The height of the density around k = 1.5 fm~^ 
is related to the components of higher partial waves induced by the tensor force. The a 
particle contains the largest bump among the s-shell nuclei, but the component contained 
in the first excited 0^ state (a*) is much smaller as understood from the 3A^ + cluster 
structure. The kinetic energy density extends beyond k = 6fm~^ for all the nuclei, which is 
of course due to the short-range repulsion of the central interaction. 



IV. SUMMARY 

We have shown that the ground state energy of a nucleus as a self-bound system is a 
functional of a set of internucleon correlation functions including the kinetic energy term. 
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FIG. 5: Two-nucleon relative kinetic energy density per unit wave number. 
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Conversely the set of the internucleon correlation functions uniquely determines the nuclear 
Hamiltonian. Namely, there is a one-to-one correspondence between the nuclear Hamiltonian 
and the internucleon correlation functions. The ground state energy becomes a minimum 
for a set of the exact internucleon correlation functions. 

Using the accurate wave functions for s-shell nuclei, we have calculated the internucleon 
correlation functions for d, '^H, '^He and the ground state of ^He. To see the dependence 
of the correlation functions on nuclear structure, we have also included the first excited 
0"*" state of ^He. We used two different potentials, AV8' and GSRS, as the two-nucleon 
interaction. Both of them contain central, tensor and spin-orbit components. We have 
shown that the magnitude and the shape of each of the internucleon correlation functions is 
clearly understood from the characteristics of the underlying two-nucleon interaction. The 
correlation function for the two-nucleon relative kinetic energy also indicates the importance 
of the tensor correlation and the short-range central repulsion in its momentum dependence. 

We have discussed the asymptotic behavior of the internucleon correlation functions. For 
a large separation of the two nucleons, the correlation functions are determined by negative 
energy solutions of three-body systems interacting via Coulomb potentials. 

Studying the internucleon correlation functions for heavier nuclei will be interesting and 
important because they give us direct information on the distribution of the internucleon 
motion. We expect that the shapes of the functions do not differ drastically from those in 
the lightest nuclei, but the heights of the peaks are expected to be smaller, and the larger 
the mass number, the larger the spatial extension of the correlation function. 



APPENDIX A: CORRELATION FUNCTION FOR SPIN-ORBIT FORCE 

In this appendix we show a method of calculating the iCF for the spin-orbit force following 
the formulation of Ref. 3. We use the notation used there to be consistent with the 



formulation, so that a reader is referred to Ref. [igI] for details. Since the matrix elements 
of the spin and isospin parts are calculated using a standard method, we here focus on the 
spatial part only. 

The spatial part of the spin-orbit force has the form V{\rij\){rij x TTjj). We can express 
'^ij = fi — fj in terms of a linear combination of the relative coordinates Xk-, that is Vij = 
Ylk=i "^kXk = wx. Similarly TTjj is expressed as TTjj = Yl!k=i ^k^^k = ^tt, where tt^ = 
—id/dxk is a momentum operator conjugate to Xk- Thus the spin-orbit force V{\rij\){rij x 
TTjj) is expressed as V{\wx\){wx x ^tt). 

The spin-orbit matrix element is calculated from the following expression 

{F{L3L4)L'M'iUs, M4, A', X)\V{\WX\){WX X ^'rT)m\F(LiL2)LMiUi, U2, A, x)) 

X {g{Xse3U3 + X4eiU4]A\x)\V{\wx\){wx x'^7z)m\9iXieiUi + \2e2U2; A,x)) ,(A1) 

A,=0 

where {a x b)m stands for —^/2i[a x = — (4-\/27r/3)ia6[yi(d)Yi(6)]im. To calculate the 
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spin-orbit matrix element between the generating functions, we make use of the relation 

{g{X3e3U3 + X4,e4U4] A' , x)\V{\wx\){wx x ^7z)m\g{XieiUi + X2e2U2; A, x)) 

dr V (r) {g{XsU3 + X^e^Ui; A' , x)\5{wx - r){wx x E_TT).J\g{XieiUi + X2e2U2] A, x)) , 



(A2) 



with 



(5f(A3e3M3 + A4e4M4; A', x)\5{wx - r){wx x ^7r)|5((AieiMi + A2e2M2; A, x)) 



det5 



r X (^z + c^AB wwB v 



X exp (h)B-^v-^c{r-wB-^vy^ , (A3) 

where B = A + A', c = {wB-^w)~^, and z = Y.i=i XiCiA^B'^ with Ai = A2 = A', A3 = 
A4 = —A. When the radial form of the spin-orbit potential is scalar, i.e. is a function of 
r, we may omit cC,AB~^wwB^^v thanks to the relation 



J dr V{r){r x a) exp ^"■^^(r — a)^^ = 0, 



(A4) 



which leads to 



{g{X3e'iU'i + XieiUi\ A' ,x)\V{\wx\){wx x (^Tz)J\g{XieiUi + X2e2U2\ A,x)) 



j dr V{r){r x ^z) 



detB 

xexp I --^cr"^ + ^ pijXiXjei ■ ej + c^7iAiei ■ r j , (A5) 

\ j>i=l i=l / 

where 7^ = wB^^Ui, pij = UiB^^uj — C7j7j, and we set 

4 

(r X ^z)^ = -V2iY^ T^X^lr x e^ji^ (A6) 
«=i 

with ^ 

= eA^^-^M,. (A7) 

The symbol — > in Eq. (lASP indicates that the Af terms in the exponent, which give no 
contribution to the required matrix element, are dropped. See Ref. [l6| for details. Substi- 
tution of Eqs. (lASP and (\A6\} into Eq. (lAip yields a basic equation to obtain the spin-orbit 
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matrix element 



{F(L3U)L'M'iu3, M4, A', X)\V{\WX\){WX X ^n)m\Fi^LiL2)LM{Ul, ^2, A x)) 



(27r)^-2cV 



det5 



I drV{r)e- 



X 



LM 



X 



4 4 

-Xa[r X ec,]imexp [ ^ pijXiXjBi ■ ej + c^'jiXiBi ■ r 

\j>i=l 



i=l 



(A8) 



We now have three expressions depending on e^, two in the exponent and one in Xa[r x 
Ga]im- Using the formula (B.4) of Ref. [l6|, we can rewrite the first expression as 



exp ^ PijXiXjBi ■ Bj 



Kj>i=l 

■ E ( n + E^(P13P14P23P24;«:) 

Pij \j>i=l / K 

X Y{pi2Pl2P34P34Pl3+PuP23+P24Pl3+P23Pl'i+P24 0KK0; Kk) 

^ [^12+Pl3+Pl4(^l)^12+P23+P24(^2)] ^ [^13+P23+P34 (^3)^14+P24+P34 (^4)] ^ 
^ ^Pl2 +P13 +P14 yPl2+P2a +P24 j)^P13+P23 +P34 j)^P14 +P24+P34 



J 00 



(A9) 



The symbol ^ indicates that the angular momentum coupling must be made to its maximum 
value for each e^. See Ref. 1^ for details. The coefficients X and Y are defined in Eqs. (B.7) 
and (B.9) of Ref. 16|. With the use of Eq. (B.22) of Ref. [l^, the second term reduces to 



exp c 



^liX, 



Bi ■ r 



i=l 



j2 1 Yi tin^r^. 



5„ 



Qi \i=l 

X [[[r,,(ei)n,(e2)]^[n3(e3)n4(e4)]^'].^K^ 



5^ v/27TTC(gig2; C(g3g4; /i') C(/i/i'; 



J 00 



\qi \q2 \ 93 \q4 
^1 ^2 -^3 ^4 • 



(AlO) 



Here C is a coefficient which couples two spherical harmonics with the same argument 

C(/i/2;/3) = 



'(2/i + l)(2/2+l) 



47r(2/3 + i: 



(/1O/2OI/3O). 



(All) 



Note that i in the above summation 



Xa[r X Ba]u 



must be 1 to get a nonvanishing contribution in 
/es like a vector in r, which is simply given as 

■r[Fi(e«)Fi(r-)]i^A„. (A12) 



Eq. (lASIl because the term Xa[r x Ba]im behaves like a vector in r, which is simply given as 
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From Eqs. flAlOp and flA12D we have 



[[n.(ei)n.(e2)],[n3(e3)n4(e4)],']iri(f)J^^AfArArAr[Fi(e„)ri(f)]i^A, 
V |in,(ei)n,(e2)|jn,(e3)n,(e4)L'hn(eJ 

Im 



[[n,(eOn,(e2)]^[n3(e3)n4(e4)];.']iyi(e.) AfArA^AfA^, (A13) 



where the two Y'i(f)'s are coupled to a scalar. Note that the integration over r in Eq. (lASP 
gives 4:71. It is convenient to express Eq. (lAlSP as 



- Im 



[[n,(eOn,(e2)],[n3(e3)n4(e4)],']in(e„)l_ArArArArA 

vv' 

x[[n..(ei)n..(e2)].[na.(e3)n4.(e4)].'],„A?-ArArA 



94a 
4 5 



where 



(A14) 



(A15) 



An expression for the coefficient will be given later. The coupling of Eqs. flAQp and 
(IA14P must lead to [[1/^^ (ei)Fi2(e2)]L [11,3(63)11,4 (64)]/^/] in order to have a nonvanishing 
contribution in Eq. flASp . This coupling gives the following factor 



2+P13+P14 



P12+P23+P24 



,62)] ^ [lpi3+P23+P34 (63)lpi4+P24+P34 (^4)] ^ 



J 00 



Im 



X0n..(ei)n2.(e2)].[n3.(e3)n4.(e4)].'] 

^Y{pi2 +Pl3 +PU Pl2 +P23 +P24 Pl3 +P23 +P34 Pu +P24 +P34 qiaq2aq3aq4a KVv'l; LL') 

X [[YLAei)YUe2)]L[YUe,)YUe,)]u],^. (A16) 
The values of pij and qi must satisfy the following equations: 

Pl2 + Pl3 + Pl4 + gia = ^1, Pl2 + P23 + P24 + ^20 = ^2, 

Pl3 + P23 + P34 + q3a = L3, Pu + P24 + P34 + q^a = ^4- (Al7) 



The operation prescribed in Eq. flASp is now easily performed. 
To sum up these results, we obtain the following formula: 

{F(LsL4)L'M'{u3, M4, A', X)\V{\WX\){WX X ^7r)„|F(iiL2)LM(Ml, ^2, A, x)) 

■(2vr) 



10 (_l)Li+L2+L+L'+l 

47i\l-- ' , {LMlm\L'M') JJ 



vi=l 



^N-2^\ 2 



det5 



X 



E^"EI n ^ '^"^y^^ (P.)"- 1 E-^(P13P14P23P24;.) 



a=l Pij \i>i=l 
X F(Pi2Pi2P34P34P13+P14P23+P24P13+P23P14+P24 0kkO; Kk) 

X E f n ^"^^y^^^" I E ^(^1^2; /i) C(g3g4; /i') 1) 



gi \i=l 

r>00 



MM 



/>oo 

X / rfr r^^+^2+^3+''*+3 e"^""' ^^(r) ^ K^{q^q2q3q4fxfJ.'; uu') 

xY{Li-qia L2-q2a ^3-g3« ^4-g4a qiaq2aq3aq4a KUu'l] LL'). 



(A18) 
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The coefficients Ka are given as follows: 

K2iqiq2q3Qil^l^'; v^') = i-T^''S^'^,'C{q2l; q22)U{n' jjll; lu)U{qiq2iyl; IJ^q22), 
K3{qiq2q3q4lJ'lJ''; vv') = {-)''' 5uf,C {q3l; q33)U{nfx' 11; Iz/')f/(g4g3'^'l; /^'gss), 
K4{qiq2q3q4lJ'IJ'']i^i^') = 5^^C{qil] qu)U{fifx' 11; W)U{q3q4i^'l; n'qu)- (A19) 

The values of u and u' are constrained by triangular relations which come from the unitary 
Racah coefficients U. Choosing \^(|iya;|) = (5(|tya;| — r)/r^ in Eq. flAlSp leads to the iCF for 
the spin-orbit force. 



APPENDIX B: ASYMPTOTICS OF INTERNUCLEON CORRELATION 

FUNCTIONS 

Here we discuss the asymptotic form of the iCF. Let Xi (instead of ri2) denote the 
relative distance vector of nucleons 1 and 2, and X2 denote the coordinate of their center 
of mass relative to the center of mass of the rest of the nucleus, which is called a nucleus 
R, containing N — 2 nucleons. When two nucleons are separated far in distance from R, all 
the nuclear forces can be neglected, and only the Coulomb interactions among them remain. 
The Hamiltonian H of the whole system thus reduces to 

h+t2 + + jT^^^Pi. + , rP2n + Hn, (Bl) 

Xi \^Xi + X2\ I - ^Xi + X2\ 

where ti and t2 are the kinetic energies, ti = —{h^/m)d'^/dx\, and t2 = —{f? I2\i)d'^ I dx\ 
with /i = UN — 2)m/N . The charge of the nucleus R is Z^e, and its internal Hamiltonian 
is denoted by Hr. Corresponding to this decomposition, the wave function \1/ for large xi 
takes the form 

XI KLs^2T^2lJRTA[^L{.Xi,X2)XS^2]lVT^2^JETn]jM.,TMT. (B2) 

LSi2Ti2IJrTr 

where ^1/ j^r^ is the normalized wave function of R with spin Jr and isospin Tr, though it may 
not be always an eigenstate of Hr. Let Ej^Tr be the energy expectation value which ^I/j^t^ 
gives, {'^ j^Tf.\HR\^ j^Ta) ■ The spin and isospin states of nucleons 1 and 2 are represented 
by X5i2 cind TjT^^. The summation labels of Eq. flB2p run over all possible angular momenta 
which satisfy the angular momentum couplings and the parity conservation as well as the 
Fermi statistics of nucleons 1 and 2. 

Using the asymptotic forms, fIBip and flB2p . in H'^ = E^^, we find that the wave function 
$L satisfies the three-body equation with only Coulomb potentials 

/ g2 Z Z \ 

(h+t2 + — + rPi. + :P2n $LAf^ = {E - Ej^Tn)^LM^. 

V ^1 \^Xi + X2\ \-^Xi + X2\ J 

(B3) 

The energy E — Ej^^t^ is however negative in contrast to a usual case 20|. Since ^ is a bound 



state wave function, ^lm^ must also be bound. A solution for large xi can be obtained by 
solving the above equation. Here we attempt to obtain an approximate solution by taking 
the leading term of the nucleon-R Coulomb potential. 
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Assuming that xi is much larger than X2 makes it possible to simplify Eq. (1B3P to 



tl+t2 + 



Xi 



LMr 



{E - Ej^tJ^ 



LMr 



(B4) 



with 



2n) 



(B5) 

Since the coordinates Xi and X2 are now decoupled, we find a solution of the type of ^lm^ = 
[ipi{xi)(f)x{x2)]LML- The function (j)\{x2) must be bound, that is, the matrix element of t2, 
{<Pxmx\'t2\<Pxmx) y must bccome negative, the value of which is set equal to —ti^q^ /2jj,. This is 
possible only for A = 0. Thus ^lml turns out to be of form ['?/'L(35i)0o(^^2)]LA/i- Then the 
radial part of ipi = fL{xi)YL{xi) satisfies the equation 



2d 

+ 



dxf Xi dx\ 



L{L + 1) 2Kr] 2' 

n ^ 



XI 



Xi 



/l(xi)=0. 



where 



mZe^/2h^K, K= J — (^-E + Ej^t^ - 



~2^ 



(B6) 



(B7) 



This equation is the same as a scattering equation by a Coulomb potential but with negative 
energies. A solution fhixi) which decreases for large xi is given by 



(B8) 



Here W is the Whittaker function |21|, which is given using the confiuent hypergeometric 
function F by 



WaM = .^+le-iF(&-a+ 1,26 + 1; 



+ 



r(26) 



r(i + 6-a) 



^e~^F{-b - a + i, -26 + 1; z) 



(B9) 



Here F is the Gamma function. From the asymptotic form of the Whittaker function for 
large z, we have 

hixi) x^'-'^e-^^\ (BIO) 

Substituting Eqs. (IB2p and fIBSP to Eq. (|28|) gives the asymptotic form of CQ'^\r). 

The value of k depends on whether nucleons 1 and 2 are both neutrons, or protons, or a 
neutron and a proton as well as on the value of q. The case which gives a minimum value 
of K determines the behavior of iCF at large distances. 

The behavior of C{k) for small values of k is given by the Fourier transform of the right 
side of Eq. ( IBlOp , which reads 



(27r)3/2 



dxi e~ 



ik xi 



e-^^^Ydx,] 



f L + l + f] L + r] _ 3 

K^+^{k^ + K^y-^ V 2 ' 2 ' 2' 



Yr(k) 



(Bll) 
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where F is the Gauss hypergeometric series. The result is thus obtained as 



C{k) 



k 



2L 



(A;2 + «;2)2-2„ 



L + l + f] L + T) 



r 3 e\ 



H 2 



(B12) 
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